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INTRODUCTION 

This work formulates a method for the modeling of material damping characteristics in 
distributed parameter models which may be easily applied to models such as rod, plate and 
beam equations. The general linear boundary value vibration equation is modified to 
incorporate hysteresis effects represented by complex stiffness using the transfer function 
approach proposed by Golla and Hughes. The governing characteristic equations are 
decoupled through separation of variables yielding solutions similar to those of undamped 
classical theory, allowing solution of the steady state as well as transient response. 
Example problems and solutions are provided demonstrating the similarity of the solutions 
to those of the classical theories and transient responses of non-viscous systems. 


BACKGROUND 

Classical damping models assume a simple damping term proportional to velocity 
and/or strain rate such that the equations of motion decouple readily.(See for instance 
Reismann 1 and Soedel 2 ). This modeling technique is more often an equivalent viscous 
damping representation of non-viscous damping than a representation of the real structural 
damping 3 . Banks et.al. 4 have shown experimentally that hysteretic effects dominate the 
damping mechanisms in a cantilevered composite beam. However, analysis of hysteretic 
models is quite difficult. Inman 5 proposed using the Golla and Hughes viscoelastic model 
for modeling damping in composite beams, and Slater 6 extended this theory to plates. A 
simple method of representing damping mechanisms over a wide frequency range which 
can be included in the general linear boundary value vibration equation to arrive at a closed 
form solution for the damped response is presented in the following. 

The time hysteretic stress-strain relation is given by 


o(\t) = e(x,t)E- fg(t- s)e(\s)ds (1) 

0 


where o(x,/) is the stress, x is the vector of spacial coordinates, e(x,t) is the strain, and 
the kernel g(t-s) describes the hysteresis as developed by Christensen 7 or Banks 8 , for 
example. The general linear boundary value vibration equation incorporating viscous and 
hysteretic damping may then be written as 

w(x,t) +dw(\t)+ L 2 [l ? (x)w(x > /) + L 1 w(x,0]« f(\t) (2) 
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where w(x,f) is the displacement, d is theviscous damping coefficient, and/(x,f) is the 
forcing function. The operators E 2 and L 3 are stiffness operators; L\ represents the 
integral part of ( 1) and ( ) denotes partial differentiation with respect to time. 

In modeling damping using finite element models, Golla and Hughes 9 have shown that 
the variation of material damping and stiffness properties as a function of frequency (i.e. 
time hysteresis) can be represented by 



Pf_ + Yjl 
PS + ya + T i 


(3) 


where E* ( 5 ) is the complex modulus of elasticity, E is the static modulus of elasticity, and 
s is the Laplace domain operator. The symbols a„, p n , y„, and represent constants for 
curve fitting to match complex modulus data. Ideally, k would be infinite in order to 
perfectly match the material properties over all frequencies. However, it is up to the analyst 
to determine the value of k in order to obtain sufficient accuracy. Although this transfer 
function has been developed for finite element models, it is easily extendable to distributed 
parameter theory. 

Using (3) to model the time hysteretic stress-strain effects of (1) gives 


o{x,s) = E’(s)e(x,s) = E 1 



ps + 1£ 

P/ + yfi + v, 


e(s) 


(4) 


Incorporating the GHM model into the Laplace domain representation of (2) results in 


s 2 VV(x, 5 ) - sw( x, 0) - w, (x, 0) + sdW{ x, s) 
-d\\{\, 0) + I^[E' (s)N(x)W(x, j)] - F(x,s) 


where the operator L\ has been replaced by 




P/ j Y£ 

P / + YJ + Vi 


( 6 ) 


and the operator L 3 has been factored into 

E, = EN(x) (7) 

to yield a linear partial differential equation. The purpose of the factorization of (7) is for 
the case of the beam or rod of varying cross section. For example, for a beam with varying 
cross section. 


N{x)- 


/(•*) 

pA d 1 x 


(8) 


where x is the distance along the beam, /(*) is the moment of area, p is the material 
density, and A it the cross sectional area. Note that when the force or displacement is 
represented by capital letters it denotes the Laplace transformed variable. 
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ANALYSIS 


Consider the solution of (5) by separation of variables/modal expansion with associated 
permissible boundary conditions. The solutions are then in the form 

00 CD 

W(x,5) = 2tv„(x,5)=Ja„(i)iI>„(x) (9) 

n - ! n-1 

Substituting (9) into (5) yields 

2[V + sd + £*(-s)L 2 W(x))A i (.s)]<I> b (x) - (s + rf)w(x,0) - w(x,0) = F(x, s) (10) 

n-1 J 

If <J>„(x) are chosen to be the eigenfunctions of the operator LjjV(x) such that the 
boundary conditions are satisfied and 

4^(X)0 =AcP n (11) 

then (10) becomes 

00 

2) IV + sd + E\s) A„ )A„ ($)]<!>„( x) ~{s + d)w(x, 0) - w(x,0) = F(x,s) (12) 

n-1 L J 

where A n are the eigenvalues of the associated with the eigenfunctions O n (x). Note that 
(11) is the usual modal assumption for the undamped equations of motion. 

The functions are assumed to be orthonormal such that 

j*0„(x)0 m (x)dQ= (13) 

Q 

where Q is the spacial domain of the equation, i.e. length for a beam and area for a plate. 

Multiplying (12) by O m (x) , integrating over the region Q and applying (13) yields the 
n temporal equations 

[V + ** + tf^A.K (*)-(* +rfk(0)- «„(0)]= f H (s) (14) 

where 

f n (s) =jF(x,s)<t> n (x)dQ (15) 

and 


a„(0)“ f w(\0)<P n (x)dQ 

Q 


a n (0)=Jw(\0)<P n ( x)dQ 

a 


(16) 
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The solution of the temporal equation is then 


* 



)+ 1 


where (3) has been substituted into the solution and d, = fis 2 + y,s + rj, . 


(17) 


EXAMPLE 

Consider a problem involving decay from initial conditions of a simply supported 1 cm 
thick 1.5m x 2m plate given an assumed complex modulus. A complex modulus 
represented by a single expansion term has been arbitrarily chosen in order to obtain a 
simple solution. The parameters are as labeled. 

x A 


1.5 


Figure 1 

p = 2700-^f-, £ = 6.89 x 10 10 Pa, v =.34, h =.01m 

971 

The viscoelastic properties are a =.6,(5- l,y - 1.5x 10 4 ,p = 1.5x 10 6 . 

xjz yjz 

The initial conditions are given as w(x, y,0) =.01sin(— )sin(— ) , which is the first 

mode shape, with an amplitude of .01. The following time derivatives are also set: 
w(^;y,0)- vvfjr.y.O) « w(jr,y,0) = 0. The need for this will become apparent. The 
equation of motion for a plate is 
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Z. 

a{t) =[ 5 . 6 x 10 4 5.6 x 10 4 3.7x10-2 0 ] ^ (E6) 

Z«. 

where the states z are related to the displacement by equation (E6) and the initial conditions 
are a(0)=.01 and a(0) - a(0) - a(0) =0. 

From eqn. (E6), its time derivatives, and the initial conditions, 

5.6 xlO 4 5.6 xlO 2 3.7 xlO' 2 0 lkl oi 

0 5.6 xlO 4 5.6 xlO 2 3.7 xlO' 2 Z 2 0 

-6.3 xlO 8 -1.0 xlO 7 -67 x 10 2 0 Z 3 “ 0 

0 -6.3 x 10 8 -1.0 xlO 7 -6.7 x 10 2 7 . 0 

the solution of which is 
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(E8) 




4 . 8 x 10' 7 

*2 


-16 xlO ' 6 

Z3 


-4.3x10-' 

Z 4 

r-0 

6.4 xlO 3 


The results of a digital simulation of the plate response are shown in figure 2. The decay 
envelope for an equivalent viscous damped plate = r\/2 where rj is the loss factor of the 
material at 124. 14 rad/sec) is shown for comparison. 



Figure 2 


DISCUSSION 

A method has been presented solving viscoelastic plate problems which accounts for 
frequency dependent modulus of elasticity. The solution yields the same orthogonal 
eigenfunctions/modes as classical plate theory while producing decaying temporal functions 
representing viscoelastic effects. This method provides a simple approximation for 
modelling viscoelastic plates. The solution reduces to the classical Sophie Germain 
solution w hen the hysteresis terms are dropped, as can be seen in equation ( 17) when a n is 
set equal to zero. 


552 




ACKNOWLEDGEMENTS 


This work was supported by NASA Langley Research Center, Grant No. NGT 50541, 
under the direction of Dr. Garnett Homer. 


REFERENCES 

1) Reismann, H., Elastic Plates, Theory and Application , John Wiley & Sons, New 
York, 1988. 

2) Soedel, W., Vibration of Shells and Plates, Marcel Dekker, Inc., New York, 1981. 

3) Cudney, H.H. and Inman, D.J., “Determining Damping Mechanisms in a Composite 
Beam by Experimental Modal Analysis,” International Journal of Analytical and 
Experimental Modal Analysis , Vol. 4, No. 4, Oct. 1989, pp. 138-143. 

4) Banks, H.T., Fabiano, R.H., Wang, Y., Inman, D.J., and Cudney, H.H., “Spatial 
Versus Time Hysteresis in Damping Mechanisms,” Proceedings of the 27th IEEE CDC, 
1988, pp. 1674-1677. 

5) Inman, D.J., “Vibration Analysis of Viscoelastic Beams by Separation of Variables 
and Model Analysis,” Mechanics Research Communications , Vol. 16, Apr. 1989, pp. 
213-218. 

6) Slater, J.C. and Inman, D.J., “Transfer Function Modeling of Damping Mechanisms 
in Viscoelastic Plates,” 32nd Structures, Structural Dynamics, and Materials Conference, 
Apr. 1991, pp. 2381-2383. 

7) Christensen, R.M., Theory of Viscoelasticity: An Introduction, Academic Press, New 
York, 1982. 

8) Banks, H.T. and Inman, D.J., “On Damping Mechanisms in Beams,” No. #89-3, 
University of Southern California, Sep. 1989. 

9) Golla, D.F. and Hughes, P.C., “Dynamics of Viscoelastic Structures — a Time 
Domain Finite Element Formulation,” ASME Journal of Applied Mechanics , Vol. 52, Dec. 
1985, pp. 897-906. 


553 




